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Abstract In comparison to 3d or 4/ metals, magnetism in actinides remains poorly 
understood due to experimental complications and the exotic behavior of the 5/ 
states. In particular, plutonium metal is most especially vexing. Over the last five 
decades theories proposed the presence of either ordered or disordered local mo- 
ments at low temperatures. However, experiments such as magnetic susceptibility, 
electrical resistivity, nuclear magnetic resonance, specific heat, and elastic and in- 
elastic neutron scattering show no evidence for ordered or disordered magnetic mo- 
ments in any of the six phases of plutonium. Beyond plutonium, the magnetic struc- 
ture of other actinides is an active area of research given that temperature, pressure, 
and chemistry can quickly alter the magnetic structure of the 5/ states. For instance, 
curium metal has an exceedingly large spin polarization that results in a total mo- 
ment of ^8 /is/atom, which influences the phase stability of the metal. Insight in 
the actinide ground state can be obtained from core-level x-ray absorption spec- 
troscopy (XAS) and electron energy-loss spectroscopy (EELS). A sum rule relates 
the branching ratio of the core-level spectra measured by XAS or EELS to the ex- 
pectation value of the angular part of the spin-orbit interaction. 



1 Introduction 

Actinide elements such as uranium and plutonium often conjure visions of radiation- 
spawned monsters who wreak havoc on unsuspecting bystanders. While this makes 
for spectacular pulp fiction, it does not relay the fact that actinide materials are 
the cornerstone of nuclear energy and are scientifically fascinating. At present, nu- 
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clear power offers the only viable carbon-free energy production that can meet the 
growing demands of the world that will stress the power grid over the next several 
decades. Scientifically, the 5/ electron states of actinide materials create a myr- 
iad of exciting physical behaviors, such as superconductivity, itinerant magnetism, 
electron correlations, and heavy fermions [43]. Therefore, these heavyweights of the 
Periodic Table are relevant both technologically and scientifically, posed to play a 
key role in the resurgence of nuclear power for clean and robust energy and expand 
our understanding of condensed-matter physics. 

The actinide series consists of fourteen elements and, accordingly, the 5 / states 
can contain up to the same number of electrons. Adhering to Hund's rule of maxi- 
mizing total spin, this means the /-shell accepts seven electrons with the same spin 
direction before taking electrons with opposite spin. Therefore, magnetism has the 
potential to be strong near the middle of the series. Indeed, this is observed in the 
rare-earth series where the middle elements such as Eu and Gd exhibit very strong 
magnetic moments due to spin polarization of the 4/ states. The strong magnetic 
moments observed in the rare-earth series are a direct consequence of the 4 / states 
being localized, i.e., reduced wavefunction overlap with non- or weakly-bonding 4/ 
electrons. The opposite is true for the light actinides, such as Th, Pa, U, Np, and a- 
Pu, because the 5 / states are delocalized in the metallic form. This changes in the 
latter part of the actinide series where the 5/ states become localized. The actinide 
series, therefore, does not simply adhere to Hund's rule. Rather, there is a competi- 
tion between the exchange interaction that drives maximization of the total spin and 
the spin-orbit interaction of the j = 5/2 and 7/2 angular momentum levels of the 5/ 
states [44-46, 80]. This competition produces interesting magnetic behavior in the 
actinides when the states localize in heavier elements and when the light actinide 
elements bond with other elements. 

In this Chapter, we discuss in general terms the magnetism of the actinide ele- 
ments, in particular Pu. A comparison will be draw with the 3d, Ad, and 5d metal 
series as well as the Af metals. Because electron localization is a requirement to 
obtain local magnetic moments, we will pay specific attention to this in both metals 
and some compounds. For instance, the transitions from a- to S-Pu and from a- to 
y-Ce both exhibit a large change in volume, indicative for a transition from itiner- 
ant to localized behavior. However, unlike Ce metal that shows a magnetic moment 
in the large-volume y-Ce phase with localized Af states, both a- to 5-Pu show no 
evidence of a magnetic moment. Surprises like this, litter the actinide series, often 
creating controversies that are fodder for fruitful and entertaining arguments. We 
will specifically discuss subjects such as 

• Localization of the 5 / states across series, 

• Crystal phase transition and comparison to Cerium, 

• Magnetic properties, 

• Consequences of angular-momentum coupling of the 5 / states: LS, jj, and inter- 
mediate coupling, 

• Role of spin-orbit and electrostatic interactions, 

• Electron energy-loss spectroscopy (EELS) in a transmission electron microscope 
(TEM), 
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Fig. 1 Wigner-Seitz radius of 
each metal as a function of 
atomic number Z for the 5d, 
4f, and 5/ metal series (from 
Ref. [10]). The upper-left in- 
sets schematically illustrate 
the localized and delocalized 
5/ states between adjacent ac- 
tinide atoms (from Ref. [2]). 
Note the parabolic shape of 
the 3d series, the almost con- 
stant values of the trivalent 
rare earths (for clarity the 
divalent Eu and Yb are omit- 
ted), and the unusual behavior 
of the 5/ elements, with min- 
imum volume near Pu, and 
strong increase between Pu 
and Am. 

• ^4.5 branching ratio, 

• Spin-orbit sum rule and its validity for the 5/ states of actinides, 

• Spin-orbit interaction per hole for Th, U, Np, Pu, Am, and Cm, 

• Electron population of the /j ^ an d fi /2 levels, 

• Ramifications of these results for the magnetic structure of actinides. 

We point the reader to other review articles, such as Challenges in plutonium 
science, Vol. I and II in the Los Alamos Science Series [1], particularly the article 
entitled Plutonium Condensed-Matter Physics. A Survey of Theory and Experiment 
by Boring and Smith [10], Absence of Magnetic Moments in Plutonium by Lashley 
et al. [33], Handbook on the Physics and Chemistry of the Actinides by Freeman 
and Lander [18], and Nature of the 5f States in Actinide Metals by Moore and van 
der Laan [43]. These reviews cover the electronic, magnetic, and crystal structure 
of the actinide materials. Here, we will focus on the magnetic structure of actinide 
metals and materials, incorporating the most current experimental and theoretical 
results. 



2 Volume Change Across the Actinide Series 

Figure[T]shows a plot of the metallic radii, which are directly proportional to volume 
of the actinide 5/ metals, in comparison to those of 5d and 4f metals. The curve 
for the 5d transition-metal series is parabolic because in the first half of the series 
each additional electron contributes to the cohesive energy of the solid, resulting in a 
decrease of volume until the shell is approximately half full. In the second half of the 
d series there is an increase of volume due to filling of the antibonding states. This 
behavior is indicative of states that are delocalized and actively bonding. The early 
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actinides (Th to Np) exhibit a contraction with increasing electron count similar to 
the 5d series, and this in conjunction with the absence of magnetic order indicates 
that the 5/ electrons of the early actinide elements are itinerant (delocalized). 

In the vicinity of Pu there is an abrupt jump in the volume and the elements Am 
to Cf exhibit no change in volume with increased 5/ electron count. This behavior 
can be explained by examining the 4/ rare-earth series in Fig. [TJ Here, there is no 
change in metallic radii across the series because the 4/ states are localized and 
non-bonding, leaving cohesion to the (spd) 3 electrons that do not change count 
across the series. Comparison with the 4/ elements, together with the presence of 
ordered magnetism in Cm and elements beyond, indicate localized behavior of the 
5/ electrons. 

Our discussion of the two types of bonding behaviors in the actinide series em- 
phases the unique position of plutonium. The volume change between a-Pu and Am 
is ^50%, a staggering change in volume between two neighboring elements in the 
Periodic Table considering that the only apparent change is to add one electron in the 
5/ shell. Indeed, unlike the lanthanide elements Eu and Yb, which are both divalent 
in the normally trivalent lanthanide series, there is no indication of a straightforward 
valence change between Pu and Am. This means that the simple addition of another 
5/ electron is not the pressing issue, but rather the change from itinerant to local- 
ized behavior. What we will find is that the bulk of the transition from itinerant to 
localized 5/ states occurs within Pu itself! Compared to the ground state a-phase, 
the high-temperature 5-Pu phase is ^25% larger in volume. Thus, plutonium is the 
nexus between itinerant- and localized-electron behavior. 



2.1 Photoemission spectroscopy's two cents 

Clear support for increased localization across actinide series is provided by 4f core 
and 5/ valence-band photoemission of actinide metals. When examining the pho- 
toemission spectra for the elements traversing the series, the 4/ photoemission of a- 
Th [49], a-U [49], and a-Np [50] is dominated by well-screened peaks, character- 
istic for itinerant behavior. Between a-Pu and 8-Pu poorly-screened peaks that are 
due to more localized 5/ states increase in intensity compared to the well-screened 
peak [7]. a-Am shows almost uniquely unscreened peaks, which are characteris- 
tic for localized behavior [51]. In a similar fashion, valence -band photoemission of 
actinide metals shows an increase in structure from Np to Am, where growing struc- 
ture indicates localization of the 5/ states [9,21,50,51]. The increase of structure in 
4/ core and 5/ valence-band photoemission from Np to Am shows that the entire 
transition from itinerant to localized 5/ states occurs over several elments, however, 
the a-Pu and 5-Pu spectra clearly reveal that the brunt of the change occurs within 
the metal's phases. 

The most important insight from the previous discussion is that the 5/ states in 
Th to Pu are slightly overlapping between neighboring actinide atoms. They there- 
fore occupy a very narrow energy band with very high density of states near the 
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Fig. 2 Length change of Pu 
as a function of temperature, 
including the liquid phase. 
The crystal structure of all 
six solid allotropic phases 
of the metal is given in the 
lower-right side of the figure 
(after Ref. [24]). The structure 
changes from low-symmetry 
monoclinic (a-phase) to 
high-symmetry fee (5-phase), 
which occurs with exceeding 
large volume changes over a 
short temperature span. Note 
also the 5-phase contracts as 
it is heated, and plutonium 
contracts as it melts. 
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Fermi energy. As the number of 5/ electrons populating the 5/ band increases, 
the specific properties of the band begin to dominate the bonding properties of 
the metal. At plutonium, the 5/ states rapidly change towards localization and this 
abrupt change produces a metal with more solid allotropic phases than any other 
element in the Periodic Table. 



3 The Six Crystal Allotropes of Pu Metal 

Plutonium is often described as "A physicist's dream but an engineer's nightmare" 
because the metal has an anomalously low melting point, a high number of al- 
lotropes, negative coefficient of thermal expansion in some phases, and the unusual 
preference for low-symmetry crystal structures. At ambient pressure Pu exhibits six 
solid allotropic crystal structures (a, j3, y, 8, 8\ e) between room temperature and 
melting, as shown in Fig [2] The phase stability of Pu as a function of tempera- 
ture clearly illustrates the change in volume that occurs between the low-symmetry 
monoclinic structure of a-Pu with 16 atoms in the unit cell and the high-symmetry 
structure of 5-Pu. Indeed, despite a cubic structure, the fee 5-phase is perplexing 
given its unusually low density and negative thermal expansion coefficient, i.e., it 
contracts when heated. It is, however, most intriguing that the a, fi, and 7 phases 
of Pu are low symmetry monoclinic or orthorhombic structures, a geometry that is 
common to minerals but foreign to metals and alloys which are usually cubic or 
hexagonal. An immediate question is why do the light actinides have low symmetry 
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crystal structures if most all other metals have high symmetry atomic geometries? 
This is quite important given that the high symmetry structures of most metals, 
such as fee copper and aluminium, exhibit ductility by readily allowing disloca- 
tion motion through 12 slip systems. It is the ductile behavior of metals that makes 
them prized for many technological applications. Accordingly, the low symmetry 
and concomitant brittle nature of Pu and other light actinides is unwanted and we 
must understand its physical origins. 



3.1 Lowering the Electronic Energy Through a Peierls-like 
Distortion 

Protactinium, uranium, neptunium, and plutonium all exhibit low-symmetry ground- 
state structures rather than the high-symmetry structures usually found in metals. In 
the early days of actinide research, this low symmetry was attributed to directional 
or covalent-like bonding resulting from the angular characteristics of / electrons. 
However, this so-called chemist view that the pointy /-electron atomic orbitals are 
forming directional bonds, hence no closed-packed structure of the metals, falls 
short of explaining the abundance of crystal structures observed of the actinide ele- 
ments. Accordingly, physicists looked for a better answer and found it in the width 
of the band. Soderlind et al. [65] showed that as the bandwidth of a bonding state 
is narrowed the crystal structure of the metal will distort through a Peierls-like dis- 
tortion to a low-symmetry geometry. The original Peierls-distortion model was for- 
mulated in a one-dimensional lattice: A row of equidistant atoms can lower its total 
energy by forming pairs. The lower periodicity causes the degenerate energy lev- 
els to split into two bands with lower and higher energies. The electrons occupy 
the lower levels, so that the distortion increases the bonding and reduces the total 
energy. In one-dimensional systems, the distortion opens up an energy gap at the 
Fermi level making the system an insulator. However, in the higher dimensional 
systems, the material remains a metal after the distortion because other Bloch states 
are filling the gap. This mechanism is very effective if there are many degenerate 
levels near the Fermi level — that is, if the energy bands are narrow with a large 
density-of-states [10]. Symmetry breaking of a crystal structure through a Peierls- 
like distortion due to narrow bandwidths (and concomitant high density-of-states 
near the Fermi level) is true for s, p, d, and / states, meaning that even aluminium 
can be expanded enough that a low-symmetry body-centered tetragonal structure is 
observed [65]. Of course for metals with s, p, and d bands actively bonding, huge 
theoretical tensions must be applied to narrow the bands enough to lower the crystal 
structure symmetry. However, the 5 / bands of Pa, U, Np, and Pu that are bonding 
are narrow enough to yield low symmetry structures at ambient conditions. Thus, 
the narrow width of the 5 / band causes the low-symmetry structures of the ground 
state, without the need to invoke an explanation involving directional bonding. 

Our discussion of Pu and energy-lowering lattice distortions means that the metal 
is directly on the border between between itinerant and localized 5 / electron behav- 
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ior. The metal shows a monoclinic ground state and an fee high-temperature phase, 
resulting in 5/ states that are balanced between a Peierls distortion to a low sym- 
metry structure due to stronger bonding of the / states and a high-symmetry crystal 
structure with strong electron correlations with weaker bonding of the / states. The 
fact that <5-Pu exhibits a high symmetry fee structure even though it has active, albeit 
weaker, 5 / electron bonding is somewhat perplexing. Surely the degree of bond- 
ing is less than in a-Pu given that the ground-state phase is monoclinic. However, 
the exceedingly high electronic specific heat (Sommerfield coefficient) of 35-64 mJ 
K~ 2 mol -1 (Refs. [29,34,85]) tells us that there is appreciable electron weight at 
the Fermi level in 5-Pu, and this could only be achieved with the 5 / states. In order 
to understand this apparent contradiction, we can turn our attention to a rare-earth 
metal that has 4 / states which can be changed from bonding to non-bonding through 
temperature or pressure, all the while retaining its high symmetry fee structure. 



3.2 Comparison with Cerium 

A material with similar properties as plutonium, and of comparable complexity, is 
cerium, the first element in the row of the lanthanides which is known for the am- 
bivalent character of its Af states. It exhibits multiple crystallographic phases that 
are strong functions of temperature, pressure, and chemistry [31]. At ambient pres- 
sure, Ce metal exhibits four allotropic phases between absolute zero and its melting 
temperature at 1071 K, namely a (fee), j3 (dhep), y (fee), and 8 (bcc). Magnetic, 
nonmagnetic and superconducting behavior is also observed. There is broad con- 
sensus that the unusual properties of Ce and its compounds originate essentially 
from the interplay of strong electronic correlations between the Ce 4/ electrons 
and hybridization between Af- and conduction-electron states. Phenomena such 
as intermediate-valence, which points to non-integer occupation of the Af shell, 
or heavy-fermion behavior characterized by an extremely large contribution of the 
electronic specific heat, are prominent observations on such systems. But perhaps 
one of the most intriguing phenomena is the y-a phase transition in Ce metal [31]. 
It occurs under pressure at room temperature or, at ambient pressure, on cooling to 
low temperature. Complexity of the underlying mechanisms is reflected in the col- 
lapse of the atomic volume by <~17% preserving the fee lattice symmetry, and the 
loss of the magnetic moment (y-Ce is paramagnetic and a-Ce non-magnetic). The 
y-a phase boundary ends at a critical point where the two phases become indistin- 
guishable. 

The anomalous behavior of Ce at the y-a phase transition is commonly de- 
scribed by mechanisms involving either a Mott transition [30] or Kondo hybridiza- 
tion [3,35]. However, recent first-principle calculations provide evidence for a com- 
bination of effects, involving strongly correlated Af electrons in both phases and 
screening as important ingredients [16, 39, 87]. The transition to the a-phase is 
accompanied by a sizable increase of the hybridization between Af- and valence- 
electron states. It leads to a partial derealization of the Af states, i.e., the ground 
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state occupation number, n /, that is almost equal to one in the y-phase is reduced in 
the a-phase. This has been directly derived from recent measurements of resonant 
inelastic x-ray scattering on elemental Ce under pressure [55]. The experiments per- 
mit to conclude that the changes of the electronic structure at the phase transition 
mainly result from band formation of Af electrons that concurs with reduced elec- 
tron correlation and increased Kondo screening. The experiments also highlight the 
importance of double occupancy of the 4 / states in the ground state for understand- 
ing the effects of electron correlation in this element, as stated in Ref. [39]. Con- 
sidering lattice dynamics, phonon density of states of a and y-phases in Ceo.9Tho.i 
by Manley et al. show only a very small difference in the vibrational entropy of 
the two phases, where most of the transition entropy can be accounted for with the 
crystal field and changes in the ground-state spin fluctuations [36]. 

In compounds, Ce shows y-a-phase-like character depending on the Af- valence 
band hybridization strength. Interestingly, in a-like compounds, such as CeFe2 or 
CeCos, the Ce atoms carry an ordered magnetic Af moment [11] in contradistinc- 
tion to elemental a-phase Ce metal [31]. An important concept in the theoretical 
description of these compounds is the hybridization between the 3d states of the 
transition metal and the delocalized Ce Af states. Ce readily reacts with hydrogen. 
As with other light rare earth metals, it forms a cubic dihydride CeH2 with hydro- 
gen atoms on tetrahedrally coordinated sites (CaF2-type structure) and dissolves, 
in a single phase, further hydrogen on octahedral sites up to the cubic trihydride 
CeH3 (BiF3-type structure) [73]. As the composition of the trihydride is approached 
a metal-to-insulator transition occurs that is reversible. It was recognized that thin 
films of such rare-earth hydrides allow to rapidly switch between the contrasting op- 
tical properties of the dihydride and trihydride phases, making them technologically 
interesting as optical switches [28]. Model calculations of the rare earth hydrides 
in an ionic picture have shown the importance of electron correlations for this ef- 
fect [17,52,53]. Strong Coulomb interactions between the electrons on hydrogen 
sites were shown to be responsible for opening up a gap of ^2 eV between the va- 
lence bands derived from rare earth-hydrogen and hydrogen-hydrogen hybridization 
and a set of bands predominantly of rare earth-metal d character. X-ray absorption 
spectra have shown that the electronic configuration of the Ce Af states in the Ce 
hydrides is similar as in the y-phase of the pure metal [5]. Below 7 K, CeH^ orders 
magnetically in a complex phase diagram [8,73]. 



3.3 Stabilized 8-Plutonium 

Now back to the point of why we introduced Ce. a-Ce has actively bonding / states, 
but the crystal structure is fee, as observed for 5-Pu. For a-Ce, this means there is 
just enough / bonding to influence the physical properties such as electronic spe- 
cific heat, but not enough to reduce the crystal structure. Of course this can easily be 
changed with external influences. When pressurized to about 4 GPa, the structure 
changes to orthorhombic Cmcm (the same as a-U) meaning the Af states become 
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decisive enough in bonding so as to break the cubic symmetry of the crystal. Con- 
sidering 5-Pu, a similar situation is occurring where there is enough / bonding to 
influence the electronic specific heat, but not enough to break cubic symmetry. It 
is possible there is a slight distortion from a perfect cubic symmetry in 5-Pu, but 
if true this would be small [42,74]. Once 5-Pu is pressurized to about 0.1 GPa, the 
structure looses cubic symmetry. In a reverse manner, adding a few % of Ga, Al, Ce, 
Am, or other tetravalent elements stabilizes 5-Pu at room temperature and below. 

The addition of impurity atoms often destroys the coherence of the / band. With- 
out its narrow / band actively bonding, plutonium can no longer reduce its energy 
by lowering its symmetry through Peierls-like distortion to the a -phase; it therefore 
remains in the 5 -phase. Another view is that plutonium atoms relax and move to- 
ward the smaller non-/ atoms, thereby reducing the /-/ interactions that stabilize 
the a-phase. Either way, we see that addition of many different elements stabilizes 
5-Pu, such as Al, Ce, Am, and, most often utilized, Ga. Indeed, we find 5 -phase sta- 
bilization occurs even for defects, where it has been found that excessive plastic de- 
formation, and, concomitantly the introduction of dislocations, stabilizes 5-Pu [20]. 
Once more we see how external influences, such as temperature and chemistry as 
well as the pressure or defects can influence the phase stability of plutonium. What 
is more, we now have a way to overcome the crippling technological disadvantage 
of monoclinic Pu that is brittle and difficult to machine. 



4 Revised View of the Periodic Table 

Plutonium is not alone in juggling its electrons between bonding and localized 
states. Figure [5] shows a rearranged Periodic Table containing the d and / elec- 
tron series, where the / series are placed on top and spaced slightly tighter. The 
localization increases going from bottom to top. Furthermore, along each series the 
localization increases from left to right. The elements on the white diagonal stripe 
form the rough dividing line between localized (local moment magnetism) and itin- 
erant (superconductivity) long-range collective behavior. The elements on or near 
the transition between localized and delocalized behavior, i.e., the white diagonal 
stripe, exhibit a large number of solid allotropic phases [see Fig.|3jb)]. This is clearly 
illustrated in the lower part of Fig. [3] where the lighter shades indicate more phases, 
and the light band mirrors the white diagonal stripe in the rearranged Periodic Table 
above [48]. The most notable elements with numerous phases are: Pu with 6, Ce 
with 4, Mn with 4, as well as La, U, Np, Am, and Fe, each with 3. 

The elements at the diagonal white stripe in Fig.[5]have in common that the wave 
functions from different atoms are barely overlapping, and their electrons are thus 
bordering on being localized. These are the so-called correlated-electron materials 
that are characterized by the presence of a narrow conduction band, giving rise to 
spin and charge fluctuations associated with the low-energy excitations. The ex- 
otic behaviors of these materials necessitate a description by many-electron models, 
such as the Kondo, Hubbard, or Anderson models. These models can be generically 
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Fig. 3 (Color online) Rearranged Periodic Table showing the five transition metal series, 4/, 5/, 
3d, Ad, and 5d (after Refs. [62] and [10]). When cooled to the ground state, the metals in the blue 
area exhibit superconductivity while the metals in the red area exhibit magnetic moments. The 
white band running diagonally from upper left to lower right is where conduction electrons change 
from itinerant and pairing to localized and magnetic. Slight changes in temperature, pressure, or 
chemistry will move metals located on the white band to either more conductive or more magnetic 
behavior, (b) Version of (a) where the number of solid allotropic crystal structures for each metal is 
indicated by gray scale. Lighter shades indicate more phases. Notice that a band of lighter shades 
mirrors the white band in (a), showing that metals on or near the transition between magnetic and 
superconductive behavior exhibit numerous crystal phases. 



classed as two-electron "impurity" models, which introduce interactions between 
pairs of electrons, one localized on an impurity atom and one in a conduction band. 
The Anderson Hamiltonian contains both a repulsive Coulomb term and an electron 
hopping term. The former term kepdf the 5/ electrons localized, and the latter term 
leads to a partial localization of the conduction electrons. [10] 



5 Actinide Magnetism 

The relation between /-electron localization and the formation of magnetic mo- 
ments, as illustrated in Fig. [3|a), leads to the following conclusions. The light ac- 
tinides (Th, U, Np) are delocalized and show no magnetic ordering. Many of these 
metals are even superconducting at low temperatures. On the other hand, the local- 
ized heavy actinides (from Cm onwards) show magnetic ordering. Note that ameri- 
cium is 5f 6 with a ground state 7 = and therefore has no magnetic ordering. Thus, 
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Fig. 4 Molar susceptibility 
of Mn and various forms of 
Pu plotted as a function of 
the melting point [r me it(Mn) 
=1519 K and (r mdt (Pu) = 913 
K] . The susceptibility of Pu is 
characteristic of metals with 
relatively strong paramag- 
netism. Manganese is given 
for comparison. The sequence 
for Mn is a (complex cubic), 
j8 (complex cubic), y (fee), 
and 5 (bcc). Mn orders an- 
tiferromagnetically at 95 K 
(T/T me i t = 0.063) in a com- 
plex tetragonal structure that 
is a slight distrotion of the 
a-phase. From Ref. [33]. 
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magnetism — or lack thereof — of most of the actinides is clear. The situation regard- 
ing the magnetism of Pu, however, remains murky. 



5.1 Experimental Absence of Magnetic Moments in Plutonium 

Plutonium shows a large experimentally observed magnetic susceptibility, hints of 
heavy-fermion behavior, and an anomalous temperature dependence of the elec- 
trical resistivity. At low temperatures the resistivity increases, which is character- 
istic for a Kondo lattice. Judging by its small volume the monoclinic a-phase is 
clearly delocalized and therefore non-magnetic. The fee 5-phase is more localized. 
So, could this phase be magnetic? In the past, ab-initio electronic structure calcu- 
lations that included spin and orbital polarization successfully reproduced the large 
volume change by localizing the 5/ electrons in 5-Pu [4,32,57,64]. One particular 
set of calculations by Soderlind and Sadigh [67] achieved the appropriately spaced 
energies and atomic volumes for five of the six allotropic phases of Pu (the high- 
temperature bcc £-phase was too high in energy) as well as an equation-of-state, 
bulk modulus, and elastic constants that are in agreement for most all of the six 
allotropic phases. However, these calculations gave rise to a new problem, namely 
they concluded that Pu contained local magnetic moments. While experimental re- 
sults in the early days were ambiguous, today it well established that all six phases 
of Pu, including 8, show no evidence of long-range magnetism. 

Magnetic measurements on Pu were reviewed in 2005 by Lashley et al. [33] 
In order to understand the large volume expansion between a- and 5-phases of 
plutonium a localization of the 5/ states is needed. With strong reasons for taking 
the number of 5/ electrons in Pu as five, this gives a Kramer's ground state with 
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the prediction of magnetic ordering for the 5-phase, and in some cases, even for the 
a-phase. The resulting moments are predicted to be large, and even though, due to 
their antiparallel coupling there is a partial cancellation of the spin and orbital parts, 
the total magnetic moments are of the order of 1-2 /ie/atom. This remains the case 
when band structure calculations take into account orbital polarization effects, in 
which case the absolute value of the orbital magnetic moment exceeds the value of 
the spin magnetic moment. 

The measured magnetic susceptibilities of Pu in its various phases is shown in 
Fig.|4](Ref. [33]. This susceptibility is characteristic of metals with relatively strong 
paramagnetism caused by electronic band magnetism. Neither the temperature nor 
magnetic-field dependences of the measured susceptibilities in Pu provide evidence 
for disordered or ordered moments. Specific-heat measurements indicate the ab- 
sence of magnetic entropy. The combination of neutron elastic and inelastic scatter- 
ing data shows no convincing evidence for either long-range ordered or disordered 
(static or dynamic) magnetic moments. On the basis of resistivity measurements for 
both a- and 5-Pu, and review of other measurements, Arko, Brodsky, and Nellis [6] 
concluded in 1972 that there was no evidence for localized moments or magnetic 
order. On the basis of a T 2 dependence of the resistivity near T= these authors 
suggested a model involving spin fluctuations for elemental Pu, as well as for a 
number of other alloys and compounds. 

Review of the magnetic measurements by Lashley et al. [33] demonstrated be- 
yond all reasonable doubt that there is no ordered magnetism involving the 5/ elec- 
trons in Pu metal in either the a- or 5-phases down to a base temperature of ~4 K. 
The experimental evidence presented includes magnetic susceptibility, specific heat 
(with an applied field of up to 14 T), nuclear magnetic resonance, and elastic and 
inelastic neutron scattering. Previous reports of "anomalies", seen especially in the 
specific heat, can in fact be ascribed to impurities or structural effects, most prob- 
ably the occurrence of martensitic transformations of some parts of the samples to 
the a'-phase that contains Ga. The absence of any diffuse scattering in the neutron 
diffraction patterns (except at high-temperature from thermal disorder) also argues 
against any disordered local moments. Dynamical mean-field theory (DMFT) re- 
sults by Savrasov, Kotliar, and Abrahams [58] imply, that the local moments are 
"washed out" over short time scales and thus may not be observable to probes such 
as NMR and neutron inelastic scattering, depending on the probes' observational 
frequency window. Again, more effort to quantify these predictions for experiments 
on Pu would seem worthwhile. [25, 33] 



5.2 Looking to other elements for clues 

As plutonium has no local magnetic moments, a limited comparison could perhaps 
be drawn with palladium, which finds itself also on the diagonal stripe in Fig.[3ja). 
As bulk compound, palladium is an anomalous paramagnet on the verge of ferro- 
magnetism [68, 84]. The free Pd atom has a [Kr]4flf 10 closed-shell electronic con- 
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figuration. In bulk Pd no spontaneous ferromagnetic order is observed. Although 
the Fermi level is located immediately above a sharp peak in the density of states, 
at the equilibrium volume the Stoner criterion is not satisfied. The most sensitive 
factors that may enhance the density of states at the Fermi level, and cause the onset 
of ferromagnetism in Pd clusters are the electron localization associated with a re- 
duced coordination number, the expansion of the lattice, and the change of the local 
symmetry [84]. 

One might also ask whether the conduction electrons of Pu could make it an 
itinerant magnet, like iron. While there is little evidence for magnetism in pure plu- 
tonium, many plutonium compounds are indeed magnetic and tend to be itinerant 
magnets [56]. Indeed, simply dissolving hydrogen in plutonium is enough to make 
the electrons localize and to turn the system ferromagnetic. Also, a comparison of 
the light actinides with the transition metals indicates that the light actinides should 
be superconductors unless they have local moments. So, the fact that plutonium is 
not a superconductor might indicate that plutonium is an incipient, weak itinerant 
magnet and that the loss of magnetic ordering with heating plays a role in the con- 
traction of the 5 -phase. 



6 Experimental Complications of Plutonium 

Due to its toxic and radioactive nature, hands-on investigation of Pu can only be 
done at a small number of institutions around the world, making experimental re- 
sults sparse in comparison to other materials. Furthermore, experiments on Pu suf- 
fer from the problem of self-irradiation, which heats the sample, slowly destroys the 
crystal structure [41,59], grows in daughter products of the decay process, and po- 
tentially creates local magnetic moments [38]. When a plutonium nucleus in the ma- 
terial undergoes a-decay, the recoiling uranium and helium nuclei knock plutonium 
atoms from their lattice sites. Displaced plutonium atoms come to rest at interstitial 
sites and leaving lattice vacancies. Each displaced plutonium atom creates a Frenkel 
pair — consisting of a vacancy and a self-interstitial. Each decay event creates more 
than 2000 Frenkel pairs. The self-heating means that in practice it is extremely dif- 
ficult to achieve temperatures much below ^2 K, unless very small (<mg) samples 
are used, and even then there is always some doubt as to the real temperature of 
the sample. This, however, may be overcome by choosing the right isotope. The 
most common isotope is fissile 239 Pu with a 24,400 year half-life, which is made in 
nuclear reactors: 238 U + neutron — > 239 Pu. The 239-isotope produces considerable 
self heating. Alternatively, 242 Pu has a 376,000 year half-life, which drastically re- 
duces heat due to self irradiation. Finally, the creme de la creme is 244 Pu with a ~80 
million year half-life, meaning the isotope is near non-radioactive. This isotope, if 
separated in usable quantities (it is very rare), could open many new avenues for 
low-temperature experiments on Pu. 
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7 One Man's Electron Energy Loss is Another's X-Ray 
Absorption 

While the technique of x-ray absorption spectroscopy has become routine and is 
the go-to measurement for physicists, actinides and their special requirements de- 
mand a different experiment: Electron energy-loss spectroscopy (EELS) in a trans- 
mission electron microscope (TEM). Why is this? First, the TEM utilizes small 
samples, allowing one to avoid the handling of appreciable amounts of toxic and 
radioactive materials. The alternative XAS performed at a multi-user synchrotron 
radiation facility is usually less well-adapted for the delicate and secure handling of 
radioactive materials. Second, the technique is bulk sensitive due to the fact that the 
electrons traverse ^40 nm of metal, this being the appropriate thickness for qual- 
ity EELS spectra of actinide materials at the primary energy used in our TEM. A 
few nanometers of oxide do form on the surfaces of the TEM samples, but this is 
insignificant in comparison to the amount of metal sampled through transmission 
of the electron beam. Third, actinide metals near the localized-itinerant transition 
exhibit numerous crystal structures that can coexist in metastable equilibrium due 
to the small energy difference between the phases [48]. Therefore, acquiring single- 
phase samples of metals at or near this transition is uncertain, making spectroscopic 
techniques with low spatial resolution questionable. Finally, actinide metals readily 
react with hydrogen and oxygen, producing many unwanted phases in the material 
during storage or preparation for experiments. The TEM has the spatial resolution 
to image and identify secondary phases [19,27,54], ensuring examination of only 
the phase(s) of interest. A field-emission-gun TEM, such as the one used in these 
experiments, can produce an electron probe of ^5 A meaning recording spectra 
from a single phase when performing experiments is easily achieved. Quantitatively 
measuring the reflections in the electron diffraction pattern, as well as other crys- 
tallographic orientations, proves that the correct phase is examined [47, 86]. Hence 
there is no need for large single crystals! 

EELS spectra of the actinide metals have been collected at the Lawrence Liv- 
ermore National Laboratory (LLNL) using a Phillips CM300 field-emission-gun 
TEM, equipped with a Gatan imaging filter [48]. The momentum transfer in EELS 
gives rise to multipole transitions. However, by employing a small objective aper- 
ture and using a high primary energy of the electron beam (297 kV) the low-energy 
transitions are electric-dipole. The similarity between high-energy EELS spectra 
and synchrotron-radiation-based XAS has been well established for some time, and 
was recently validated for /-electron systems for the case of Ce metal [40]. 



8 Theory 



The behavior of the 5 / electrons in the actinides is governed by the interplay of the 
spin-orbit and electrostatic interactions. Here, we treat the effect of these interac- 
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tions on the electronic configuration in the different angular momentum coupling 
schemes, i.e., jj-, LS-, and intermediate coupling. This will be illustrated for the 
example of the f 2 configuration. This two-particle state is rather straightforward as 
it does not require explicit use of coefficients of fractional parentage or creation and 
annihilation operators. It is also shown how to obtain the expectation values of the 
moments in each of the different coupling schemes. For further reading on these 
topics we refer e.g. to the Hitchhiker's Guide to Multiplet Calculations [77]. 



8.1 Atomic Interactions 

For n electrons moving about a point nucleus of charge, the Hamiltonian can be 
written in the central field approximation as 



where H e i and H so are the terms for the electrostatic and spin-orbit interaction, re- 
spectively. This atomic Hamiltonian can be embedded in a solid state, e.g., using an 
Anderson impurity model [78]. Crystal field interactions are usually less prominent 
in the actinides, leading only to a small perturbation, therefore we assume a spher- 
ical potential. In any case, all interactions can be separated in angular and radial 
parts. The angular parts depend on the angular quantum numbers of the basis states 
of the electronic configuration and are independent of the radial wave functions, 
which in the calculations are taken as empirical scaled parameters (see Sec. |8.3[ ). 

The basis wave functions are assumed to be an antisymmetrized product of one- 
electron functions. In spherical symmetry, these basis states are eigenfunctions of 
the total angular momentum J and its component Mj. The states are characterized 
by quantum numbers aLS, where a is a suitable quantity for distinguishing between 
terms having the same values of the orbital and spin angular momenta, L and S. 

8.1.1 Electrostatic Interactions 

The non-relativistic Hamiltonian for the electrostatic interactions of n electrons in 
an atom with nuclear charge Ze is 



The first term describes the kinetic energy of all electrons, the second one gives 
the potential energy of all electrons in the potential of the nucleus. The third term de- 
scribes the repulsive Coulomb potential of the electron-electron interaction, where 
the expectation value can be expressed as 



H = H e i + H, 



(1) 



H e l 




(2) 
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aLS) = Y,fkF k + Y,SkG k , (3) 

Ik k 

with angular parts and gk and radial parts F k and G k . 

The radial integrals F k and G k are theoretically defined by the Slater integrals 
[12] and experimentally treated as empirically adjustable quantities to fit the ob- 
served energy levels and their intensities. The direct integrals F k represent the ac- 
tual electrostatic interaction between the two electronic densities of electrons £ and 
£'. The exchange integrals G k arise due to the quantum mechanical principle that 
fermions are indistinguishable, so that the wave function is totally antisymmetric 
with respect to permutation of the particles. Consequently, G k is not present for the 
configuration £", where the electrons are equivalent. 

The angular parts /> and can be written in terms of 3 j- and 6 j-symbols [14] 

= c-im^] (5 5 5) (o S o) {M t } » (4) 

^) = (-i) 5 [^'](ooo) 2 {^^}- (5) 

The triangle conditions require that fk(£,£) has non-zero values for k = 0,2, ...,2£, 
whereas g k {l,l') has non-zero values fork = \£-£'\,\£-£'+2\. ..,£+£', with £=££'. 
Thus the initial state f has the radial parameters F°, F 2 , F 4 , F 6 , while the final state 
d 9 f" has F 2 , F A , G\ G 3 , G 5 . 



aLS 



n-r 2 \ 



8.1.2 Spin-Orbit Interaction 

Turning to the second term of Eq. ([T]), the spin-orbit interaction for the £ shell is 
given as 

H S0 = &(r)f> Sl -, (6) 

i=l 

where 1, and s, are the one-electron orbital and spin angular momentum operators 
of the /-th electron of the £" configuration. Q(r) gives the radial part and £,1, -s, 
gives the angular part of the spin-orbit operator. For brevity, we also introduce the 
following notation 

l-s = W° = £l,.s ; , (7) 

i=i 

where the coupled tensor for the spin-orbit interaction (see e.g. Refs. [76,82] 

for details on coupled tensors). The expectation value for a one-electron state \£sj) 
with £ ^ is 



(£sj\l ■ s\£sj) = i + 1) - £(£+ 1) -s(s+ 1)] , 



(8) 



Magnetic Structure of Actinide Metals 



17 



which gives a doublet \j = £±s) with expectation values 

For the atomic many -electron state, the Hamiltonian H so commutates with J 2 and 
J z and is hence diagonal in J and independent of the magnetic quantum number Mj. 
However, it does not commute with L 2 or S 2 and can thus couple states of different 
LS quantum numbers. For a pure state \£"aLSJ) the expectation value is 

(£ n aLSJ\l-s\£"aLSJ) = hj(J+l)-L(L+ 1) -S(S+ 1)] , (10) 

where the factor £(aLS)/£ is equal to rT 1 , 0, and -(4£ + 2-n)~ 1 for n <2£ + 1, 
n = 2£+l, and n > 2^+ 1, respectively. 

It is useful to have a general expression for (1 • s) valid for any many-electron 
state in intermediate coupling, including the LS- and 7'7-coupling limits. Since the 
spin-orbit operator is diagonal in j, there are no cross terms between \ ji) and I /2). 
Therefore, we can distribute n over nu and n,- 2 , which are the electron occupation 
numbers of \ and \j2). For the configuration £" with n = rij t +nj 2 , application of 
Eq. ([8]l gives the anticipated general expression 

(f7|l-s|f7)= £ (j\l- S \j)n j = - 1 -(£+l)n ]l + 1 -£n J21 (11) 
which does not depend on the values of L, S, and J. 



8.2 LS- and j j-Coupling Schemes 

The electrostatic interaction is diagonal in the LS-coupled basis and the spin-orbit 
interaction is diagonal in the jv'-coupled basis. Switching between the LS- and jj- 
coupled basis states is done using a transformation matrix containing the recoupling 
coefficients. 



8.2.1 LS Coupling 

In the LS-coupling scheme, which gives the eigenfunctions in the limit Q(r) — > 0, 
the particles are coupled as 

[{£ a J b )L,{s a ,s b )S]J . (12) 

The Russell-Saunders notation for the L57-coupled states is 2s+l Lj. The lowest en- 
ergy term of a configuration f, i.e., the Hund's rule ground state, is that term of 
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maximum S which has the largest value of L. In addition, according to the third 
Hund's rule, J = L — S(J = L + S) for less (more) than half filled shell. 



8.2.2 jj Coupling 

In the y'y-coupling scheme, which gives the eigenfunctions in the limit {F k ,G k } — > 
0, the particles are coupled as 

[(£a,Sa)ja,(£b,Sb)jb]J ■ (13) 

In the y/'-coupled ground state, first all j = t — s levels are filled prior to the j = 
£ + s levels. For equivalent electrons, i.e., with equal values of j, the Pauli exclusion 
principle limitations on the possible values of mj prohibit some of the /-values that 
would be predicted by a vector model [14]. The result is that for any two equivalent 
electrons f 2 only the even /- values are allowed when j a = j b . 

8.2.3 Transformation Matrix 

For a given electronic configuration, we arrive in both coupling schemes at the same 
set of allowed values of the total angular momentum J = \L — S\,...,L + S. This 
means that the Hamiltonian of Eq. ([TJ is block-diagonal in /. For each /-block the 
basis states can be transformed between LS- and yj-coupling using recoupling co- 
efficients that can be expressed in terms of 9 j symbols, 

, f 4 4 l ] 

([(£ a ,ib)L,(s a ,s b )S]J\[(£ a ,s a )j a ,(£ b ,s b )j b ]J) = [L,S,j a J b ]* i s a s b S \ , (14) 

I ja jb J J 

where [x,y, ...] = (2x+ l)(2y+l) These coefficients form a matrix T = (LS \ jj). 

For a given /-block, transformation of the wave function or operator is given by 

\jj)=T\LS) and |LS)=T- 1 |07)) ) 

H (jJ) = TH^T- 1 and H<«> - T^H^T . ( ° j 



8.2.4 LS- Versus y j-Coupled Ground State: Example for f 2 

To make the angular momenta coupling more explicit, let us examine the two- 
electron case, f 2 . The possible LSJ states for this configuration are 1 Sq, x L>i, 1 G4, 
%, 3 Po.i,2, 3 F 2 .3A, 3 #4,5,6, amounting to a total of £, (2L { + 1 ) (2S,- + 1 ) = L(2/ + 1) 
= 91 Mj levels. For two equivalent / electrons, the Pauli principle allows the 
jj states (5/2,5/2)o,2,4, (5/2,7/2)^2,3,4,5,6, (7/2,7/2)o,2,4,6- This also gives a total of 
Y,i{2Ji +1) =91 Mj levels, which illustrates that the Hamiltonian in Eq. ([TJi is block 
diagonal in /. 
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In any coupling, including intermediate, the ground state has total angular mo- 
mentum J = 4. The LS-coupled Hund's rule ground state is | 3 flj) and the j j-coupled 
ground state is |(5/2, 5/2)4). Spin-orbit interaction mixes | 3 // 4 ) with | 1 G4) and 
\ 3 F 4 ). The transformation matrix for (| 3 # 4 ), l^), | 3 F 4 )) -> (|(5/2,5/2) 4 ), |(5 /2,7/2) 4 ), |(7 /2,7/2) 4 )) 



is obtained from Eq. ( 14 1 as 



L10 

3 

2 
3 



v/3 X 



V 



-2\/5 
3^2 2 



(16) 



Note that each state is (27+ 1 = 9)-fold degenerate in the absence of a magnetic field. 
Using Eq. ( [T5| , the yy-coupled ground state can be written as a sum over LS-coupled 
basis states 



|(5/2,5/2) 4 



110 



7v/3 



3 F 4 ). 



(17) 



The character is given as the square of the wave function coefficient, so that the 
/ /-coupled ground state has 74.8% \ 3 H 4 ), 22.5% 1 1 G 4 ), and 2.7% \ 3 F 4 ). This shows 
the jj -coupled ground state contains a considerable amount of low spin. 



Equation ( 15 1 also allows us to express the LS-coupled Hund's rule ground state 
as a sum over jy'-coupled basis states 



3 //4> 



l(5/2,5/2) 4 



1(5/2,7/2)4) 



1(7/2,7/2)4) . (18) 



Thus the //-coupled ground-state has 74.8% |(5/2,5/2) 4 ), 21.8% |(5/2,7/2) 4 ), 
and 3.4% |(7/2,7/2) 4 ) character. Since ( 3 H 4 | (5/2,5/2) 4 ) = 0.748, there is ~75% 
of the total electronic state can be found simultaneously in | 3 // 4 ) and |(5/2,5/2) 4 ). 
The remaining ~25% is distributed over different states depending on whether the 
spin-orbit or Coulomb interaction prevails. As we shall see, this has important con- 
sequences for the expectation values of the moments. 



8.3 Intermediate Coupling 

In order to assess the importance of intermediate coupling, we show in Table[T]the 
calculated atomic Hartree-Fock (HF) values of the radial parameters of the Slater 
integrals for representative elements among the various transition metal series. It is 
seen that the HF values of the Slater integrals F k for the different metal series are 
comparable in size. In practice an empirical scaling factor is used that depends on the 
degree of (de)localization of the valence electrons. In localized atomic systems the 
Coulomb and exchange parameters typically require a scaling to 70-80% of the HF 
value to account for interactions with configurations omitted in the calculation [72], 
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Table 1 Comparison of the radial parameters, F k (I , t ) , for the Coulomb interaction and the spin- 
orbit interaction, Q, for actinides with rare earths [72] and 2>d transition metals [79]. All values are 
in eV. The Slater integrals have been reduced to 80% of the atomic Hartree-Fock values. 







F 2 


F 4 


F 6 




25 Mn 2+ 


3d 5 


8.25 


5.13 




0.040 


64 Gd 3+ 




11.60 


7.28 


5.24 


0.197 


% Cm 3 + 


5f 


8.37 


5.46 


4.01 


0.386 



but fully itinerant systems might have to be scaled down to 1 0-20% [7 1 ] . In line with 
increasing atomic number Z, the value of the spin-orbit parameter is dramatically 
different for each of the metal series, with the largest values found for the actinides. 
While the LS-coupled Hund's rule ground state is a reasonable approximation for 
the rare earths, this will usually not hold for the actinides [15]. 

To express the wavefunctions in intermediate coupling, where {Q,F k ,G k } ^ 0, 
we can choose either the jj- or LS-basis states, but in both cases the Hamiltonian 
has off-diagonal matrix elements. The spin-orbit interaction is diagonal in jj cou- 
pling and Eq. {II} gives (1-s) = -4, -± and 3 for |(5/2,5/2) 4 ), | (5/2, 7/2) 4 ), and 



1(7/2,7/2)4), respectively, with radial part £. Transforming this diagonal Hamilto- 
nian to the LS-coupled basis (P-ffit), | 1 G4) , PF4)) using Eq 
the diagonal Coulomb interaction gives the full Hamiltonian 



15 1 and then including 



H 



(LS) 
J=4 



F( 3 ff)-3£ -J ft 



\ 

K E( l G) 
y/% £( 3 F)+|Cy 

with the energies of the Coulomb interaction obtained from Eqs. Q and Q as 



V 



(19) 
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(20) 



As is evident from the Hamiltonian in Eq. 19 the spin-orbit operator can be written 

as 

1 (21) 



L-S = LoS +-(L + S-+L-S + ), 
where the step-operators L± and mix LS states with AL = ±1 and AS = =f1. 
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Table 2 Expectation values for the spin-orbit interaction and the orbital and spin magnetic mo- 
ments for the LS-, jj-, and intermediate-coupled (IC) states of the atomic configurations f 2 and 
/ , For IC, the Slater integrals for the actinides were reduced to 80% of the Hartree-Fock values. 
Also given are the occupation numbers H5 n and n 7 n , which are related to n and (1 • s) by Eqs. |27| 
andj3§). ■ - - 
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-10 


5 





-2.86 


0.36 



8.4 Moments for f 2 

From the wavefunction coefficients we can work out the expectation value of the 
different moments. 



8.4.1 Spin-Orbit Expectation Value 



For the LS-coupled states | J // 4 ), | G4), \ 3 F4) we obtain from Eq. (19i that (1-s) = 



—3, 0, |, respectively. These results can be directly verified using Eq. (IOi, which 
is valid for pure LSJ states. 

Substituting the values of the Slater integrals F k and spin-orbit interaction £, and 
diagonalizing the Hamiltonian, we obtain the wavefunction coefficients in interme- 
diate coupling, and from those we obtain the expectation value (1-s). The result is 
shown in Tableland can be compared with the values for LS and jy'-coupled states. 
Alternatively, writing the wavefunction in yj-coupled states we can use the gen- 
An arbitrary state of the £ 2 configuration can be written 



eral expression in Eq. ( 1 1 
as 



W{t ) = cn\h,h) + cn\juh)+czL\hJ2i, 



where the wave-function coefficients fulfill the conditions n ;i = 2c 2 n 
Cj 2 + 2c\ 2 - For instance, using the expression of |3 



3 H 4 ) in Eq. 

and « 7 12 = 0.286, so that with Eq. |TT[ i we obtain (1-s) = — 2n 



(22) 

Cj 2 and tij 2 = 
18b gives n 5 / 2 = 1.714 
V --3 
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Table 3 Values of non-zero Clebsch-Gordan coefficients (LSM L M s \JMj) with \JMj) = |4, -4). 
2s+l Lj (LSM L M s \JMj) value 



3 ^4 


(5,1-5,114,-4) 


^3 




(5,l,-4,0|4,-4) 




>G 4 


(5,1,-3,-114,-4) 


v^5 


(4,0,-4,0|4,-4) 


1 


3 F 4 


(3,1,-3,-114,-4) 


1 



8.4.2 Orbital and Spin Magnetic Moments 

To obtain the magnetic moments, we assume here that the magnetic field is infinitely 
small, so that / is a good quantum number. Furthermore, we take the magnetic 
ground state at T = as Mj = —J. 

We write \LSJMj) as a sum over \LSMi,Mg) using the Clebsch-Gordan coeffi- 
cients 

(LSM L M S \J MJ ) = (-1)™^ LJ $ y (23) 

The non-zero values for f 2 (J = 4) are listed in Table [3] Only if the moments are 
stretched, i.e., L + S = J, the \LSJMj) state is a single determinant. 
The expectation values of the orbital and spin magnetic moment are 

(L z ) = (LSJMj\L z \LSJMj) = £ I {LSM L M$\JMj)\ 2 M L , (24) 
(S z ) = (LSJMj\S z \LSJMj) = £ I (LSM L M s \JMj) \ 2 M S . (25) 

M s 

The obtained diagonal elements of (L z ) and (S z ) for the LS-coupled states are given 
in Table [2] Switching to jy'-coupled basis states using Eq. ( 15 1 gives 



/ -4.57 -0.42 \ /0.57 0.42 \ 

{L z ) (ij) = -0.42 -3.8 0.54 , (S z ) ljj) = 0.42 -0.2 -0.54 . (26) 
\ 0.54 -3.43/ \ -0.54 -0.57/ 

Obviously, the trace of both (L z ) and (S z ) is conserved in different coupling 
schemes. Furthermore, the matrix (L z ) + (S z ) = (J z ) has diagonal elements equal 
to —4, while non-diagonal elements vanish. 

Results for the other /" configurations can be found in Ref. [82], which shows 
that the type of coupling has an important influence on the values of the moments. 
In Tablebl we reproduce the values of the moments for / 5 (/ = 5/2). If Pu would 
be ferro- or ferrimagnetic, then (L z ) and (S z ) could be obtained from x-ray mag- 
netic circular dichroism (XMCD). An alternative is to measure (1-s), which gives 
conclusive information about the type of angular coupling in the material. It is clear 
from from Table [5] that there is a huge variation in (1 ■ s) across the three coupling 
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cases. We will show in Sec.[lO]how to extract the expectation value of the spin-orbit 
interaction from the EELS or XAS branching ratio. 



9 Spectral Calculations 

We briefly mention how the multiplet calculations are performed. Multiplet the- 
ory [77] is ideally suited to calculate the core-level spectra for EELS or XAS at 
the A^,5 and O4.5 edges. These calculations for the transitions 5 f" — > c/ 9 5/" +1 are 
performed in the same way as for the M4 5 and Afy.s absorption edges of the rare 
earths [72], using only different radial parameters of the spin-orbit and Slater inte- 
grals (see Table[T]i. Contrary to band-structure calculations, the multiplet structure is 
calculated in intermediate coupling, which treats spin-orbit, Coulomb, and exchange 
interactions on equal footing. 

First, the initial and final state wavefunctions are calculated in intermediate cou- 
pling using the atomic Hartree-Fock method with relativistic correction [13, 14]. The 
electric-multipole transition matrix elements are calculated from the initial state to 
the final state levels of the specified configurations. At low photon energies below 
a few keV only electric-dipole transitions play a role. The electric -dipole selection 
rules from the ground state strongly limit the number of accessible final states, so 
that compared to the total manifold of final states, the allowed transitions are lo- 
cated within narrow energy regions [72], historically called white lines due to their 
appearance on photographic plates. 

Examples for the 5/° — > d 9 5f l transitions of the actinide O4.5, N45, and M45 
edges, corresponding to the 5d, Ad and 3d core levels, respectively, as well as nu- 
merical results for other f" configurations can be found in Ref. [43]. 



10 Spin-Orbit Interaction and Sum Rule Analysis 

The spin-orbit sum rule relates the angular dependent part of the spin-orbit inter- 
action to the EELS or XAS branching ratio, i.e., the intensity ratio of the core d 
spin-orbit split j'-manifolds in the f" — > d 9 f n+l transition. This is discussed in de- 
tail in this section and the validity of the sum rules is discussed in the next section. 

While the expectation value of the orbital and spin magnetic moments, (L z ) and 
(S z ), can only be measured with XMCD, the spin-orbit interaction, (I s), can be 
obtained from the branching ratio of the XAS or EELS spectra. 

For the / shell, the number of electrons, rif, and the expectation value of the 
angular part of the spin-orbit interaction, (w 110 ), is given by 



n f = « 7 /2 + «5/2 ) 

2 4 

(w 110 ) = -(l-s)=« 7/2 --n 5/2 , 



(27) 
(28) 
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Fig. 5 Calculated electric- 
dipole transitions for the A/4 5 
EELS or XAS spectra using 
many-electron atomic theory 
in intermediate coupling for 
92 U/' to/ 5 and 100 Fm/ 7 to 
/ 13 . The convolution by 2 eV 
corresponds to the intrinsic 
lifetime broadening. The 
spectrum for / 6 (not shown) 
has only a 4d 5 n peak but no 
4dy 2 peak. From Ref. [80]. 
The calculated spectra show 
very good agreement with the 
experimental A/4,5 spectra in 
Fig.[6j where n(Th) sa 1, n(U) 
fa 3, n(Np) ss 4, n(Pu) sa 5, 
n(Am) S3 6, and n(Cm) rs 7. 
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where n 7 / 2 and n.5 /2 are the electron occupation numbers of the angular-momentum 
levels j = 7/2 and j = 5/2 of the / shell. 

A sum rule [69, 70, 8 1 , 82] relates the expectation value of the angular part of the 
5/ spin-orbit interaction per hole to the branching ratio, B, of the core d to valence 
/ transitions in EELS or XAS, 

(w m ) 5 / 3\ 

4 = --(£--) , (29) 



14 - n f 2 V 5 , 

where the branching ratio for the A/4,5 edge is defined as 

B= , f 5 \ r , (30) 

with /(A/5 ) and /(A4 ) the integrated intensities of the A/5 (4t/ 5 / 2 — > 5/5 12.1/2) and A4 
(4^3/2 ~^ 5/5/2) peaks. Since the expectation value of w 110 is for the angular part of 
the spin-orbit interaction, it does not include the radial part, which is approximately 
constant for given element. Therefore, the spin-orbit sum-rule analysis is comple- 
mentary to e.g., optical spectroscopy, where energy level separations are measured. 
Thus the sum rule analysis reveals the proper angular-momentum coupling scheme 
in the material. 
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Fig. 6 Experimental N4 .5 
EELS spectra for the Ce- 
phases of Th, U, Np, Pu, Am, 
and Cm metal, normalized 
to the N5 peak height. It is 
observed that the intensity of 
the N4 (4^5/2) peak gradually 
decreases in intensity rela- 
tive to the N5 (4^3/2) peak 
going from Th to Am, then 
increases again for Cm. From 
Refs. [45,46]. The branching 
ratio of these spectra gives 
direct information about the 
expectation value of the 5/ 
spin-orbit interaction in the 
actinde metal ground state. 
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Figure[5]shows the calculated electric-dipole transitions 4d 10 5f n — > Ad 9 5 f' 1+l for 
the A^,5 EELS or XAS spectra using many-electron atomic theory in intermediate 
coupling for 92 U 5/ 1 to / 5 and 100 Fm 5/ 7 to f 13 (Ref. [82]). The line spectra are 
convoluted by 2 eV, which corresponds to the intrinsic lifetime broadening. The 
spectrum for / 6 (not shown) contains only a N5 peak while the N4 peak is zero. The 
calculated results in Fig.|5]are in very good agreement with the experimental results 
shown in Fig. [6] 



11 Validity of the Sum Rule 



The sum rule in Eq. (29i contains a correction term A, which can be calculated 
using Cowan's relativistic Hartree-Fock code [13, 14]. In the same way as the spin 
magnetic moment sum rule in x-ray magnetic circular dichroism (XMCD) [76] the 



sum rule in Eq. ( 29 1 is strictly only valid in the absence of core-valence electrostatic 
interactions, or so called jj mixing, in which case the correction term A becomes 
zero [80]. In the Periodic Table, the different transition-metal series show a differ- 
ent behavior due to diverse valence-electron interaction. For 3d transition metals, 
the application of the spin-orbit sum rule for the £2,3 branching ratio is severely 
hampered by the large 2p-3d exchange interaction that is of similar size as the 2p 
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Fig. 7 Validity of the spin- 
orbit sum rule. Theoretical 
values obtained using rela- 
tivistic atomic Hartree-Fock 
calculations show that the cor- 
rection factor, A , in Eq. [29\ is 
proportional to the ratio of the 
core-valence exchange inter- 
action and the core spin-orbit 
interaction, G l (c,£)/^ c . The 
condition for the sum rule is 
that the total angular momen- 
tum of the core level, j = c ± \ 
is a good quantum number, 
i.e., there is no mixing of 
the j = 5/2 and j = 3/2 core 
levels. This means the core- 
valence exchange interaction 
should be much smaller than 
the core spin-orbit interaction 
[G l (c,£) < £]. Numerical 
values are given in Ref. [80]. 



spin-orbit interaction [69]. The same is true for the M45 edges of the lanthanides, 
where the 3d -4/ exchange interaction is quite strong compared to the 3d core spin- 
orbit interaction [72]. On the other hand, for the L23 edges of Ad and 5a! transition 
metals, associated with the deep 2p core level, the sum rule is expected to hold quite 
well [69]. Theoretical values obtained using relativistic atomic Hartree-Fock calcu- 
lation (Cowan code) in Fig. [7] show that the correction factor, A, is proportional to 
the ratio of the core-valence exchange interaction and the core spin-orbit interac- 
tion, i.e., G l (c,£) / Thus the condition for the sum rule is that the total angular 
momentum of the core hole, j = c ± i is a good quantum number, in other words 
there should be no mixing of the j = 5/2 and j = 3/2 core states. This corresponds 
to the core-valence exchange interaction much smaller than the core spin-orbit in- 
teraction, G\c,l) < Cc [80]. 

Even in the worst case of the rare earth M45 edges, the trend in the branching 
ratio can still be used to obtain the relative population of spin-orbit split states, 
as was demonstrated for Ce systems [83]. The situation is favorable for the M4 5 
and ^4 5 edges of the actinides, given the small exchange interactions between the 
3d and 4d core levels and the 5/ valence states. This means that the EELS and 
XAS branching ratios depend almost solely on the 5/ spin-orbit expectation value 
per hole, thus affording an unambiguous probe for the 5/ spin-orbit interaction in 
actinide materials. 




La M, 



2 



Ti L. 
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12 Experimental Results for the A/4 5 Edges 

The experimental N45 (Ad — ► 5/) EELS spectra for a-phase Th, U, Np, Pu, Am, 
and Cm metals are displayed in Fig. [6] All spectra are normalized to the N5 peak 
height. Noticeable is the gradually growing separation between the N4 and N5 peaks 
from Th to Cm, in pace with the increase in 4a! spin-orbit splitting with atomic 
number. Second and more importantly, the intensity of the N4 (Ad^n) peak grad- 
ually decreases in intensity relative to the N5 (Ad^/2) peak going from Th to Am, 
then abruptly increases for Cm [45,46]. Applying the sum-rule analysis to the ex- 
perimental branching ratio gives the values of the 5/ spin-orbit interaction per hole. 

In order to visualize the spin-orbit analysis of the EELS spectra in relation to the 
results of our atomic calculations, both are shown as a plot of (w UQ )/(U-n f )~A 
versus the number of 5/ electrons in Fig. |HJa). The curves for the three theoretical 
angular-momentum coupling schemes, LS, jj, and intermediate, as calculated using 
a many-electron atomic model, are plotted as a short-dashed, long-dashed, and solid 
line, respectively. The EELS results are indicated by the blue points. Thorium metal 
falls on all three curves, due to the fact that it takes two electrons to tangle and 
with less than one 5/ electron in Th there is no difference between the coupling 
mechanisms. U falls directly on the LS-coupling curve, Np between the LS and 
intermediate curve, and Pu, Am, and Cm all fall on or near the intermediate coupling 
curve. The intermediate curve is strongly shifted towards the jj limit for Pu and 
Am, evidence of the strong preference of the 5/ electrons to occupy the j = 5/2 
level in both metals. However, for Cm there is a sudden and pronounced shift in the 
intermediate coupling curve toward the LS limit. 

The values of ny and (w 110 ) can be converted into the electron occupation num- 



bers n-j/2 and «5/2 using Eqs. (27 1 and (28 1. Note that this is just a different way of 



presenting the same data. The experimental and theoretical results are displayed in 
Fig. [8jb), where the number of electrons in the j = 5/2 and j = 7/2 levels as calcu- 
lated in intermediate coupling using the atomic model are drawn with black and red 
lines, respectively. Again, the experimental EELS results are indicated with blue 
points. Apart from the slight deviation in the lighter actinides, U and Np, which 
is caused by derealization of the 5/ states and thus indicates a departure from 
the atomic model, the EELS results are in excellent agreement with the theoretical 
curves. Figure [HJb) clearly shows that for the actinide metals up to and including 
Am, the 5/ electrons strongly prefer the f$ n level. However, this changes in a strik- 
ing manner at Cm, where not only does the electron occupation sharply increases 
for the f-j ji level, but even decreases for the f$ / 2 level. 

The physical origin of the abrupt and striking change in the values between Am 
to Cm can be understood from the angular momentum coupling. In jj coupling the 
electrons prefer to be in the f 5 1 2 level, which however can hold no more than six, so 
that the maximal energy gain in jj coupling is obtained for Am f 6 , where the / 5/ / 2 
level is full. However, for Cm f at least one electron will be relegated to the / 7 y 2 
level. The f configuration has the maximal energy stabilization due to the exchange 
interaction, with parallel spin in the half-filled shell, which can only be achieved in 
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Fig. 8 (a) Ground-state spin- 
orbit interaction per hole as a 
function of the number of 5/ 
electrons (n/), The three the- 
oretical angular-momentum 
coupling schemes are shown 
(LS, jj, and intermediate 
coupling). The dots give 
the results of the spin-orbit 
sum-rule analysis using the 
experimental W4 5 branching 
ratio of each metal in Fig. [6] 
(b) The spin-orbit interaction 
values are converted to the 
electron occupation numbers 
«5 and «7/2 of the j = 5/2 
and 7/2 levels in the 5/ state. 
Shown are the three theoreti- 
cal coupling schemes (curves 
for LS, jj, and intermediate 
coupling) together with the 
dots representing the experi- 
mental EELS results for each 
metal. [45,46] 
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LS coupling. Thus the large changes observed in the electronic and magnetic prop- 
erties of the actinides at Cm are due to the transition from optimal spin-orbit stabi- 
lization for / 6 to optimal exchange stabilization for f . In all cases, spin-orbit and 
exchange interaction compete with each other, resulting in intermediate coupling; 
however, increasing the /-count from 6 to 7 shows a clear and pronounced shift in 
the power balance in favor of the exchange interaction. The effect is in fact so strong 
that, compared to Am, not one but two electrons are transferred to the / 7 / 2 level in 
Cm [Fig. [SJb)]. Therefore, in Cm metal, the angular-momentum coupling in the 5/ 
states plays a decisive role in the formation of the magnetic moment, with Hund's 
rule coupling being the key to producing the large spin polarization that dictates the 
newly found crystal structure of Cm underpressure [23,45,66]. 



12.1 What Our Results Mean for Pu Theory 



The spin-orbit sum rule suggests an /-count near 5 for Pu, with 5.4 being a rea- 
sonable upper limit [44^-6]. Further evidence for this /-count comes from Ander- 
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son impurity calculations for Pu [75] that explain photoemission results on 1 to 9 
monolayers thin films of Pu metal [21,22]. The 4/ core-level photoemission spectra 
display a screened and unscreened peak, thereby acting as a ruler for the degree of 
localization. The results for the /-count are in agreement with recent DMFT calcu- 
lations by Shim et al. [61] and Marianetti et al. [37], as well as LDA+U calculations 
by Shick et al. [60], which explain the three-peak structure in 5/ photoemission and 
the relatively high electronic specific heat. The absence of experimentally observed 
magnetic moments in any of the six allotropic phases of Pu metal [33] is thought 
to be due to Kondo screening and a non-single Slater-determinant ground state in 
these respective models. The lack of local moments 5-Pu has also been recently ex- 
plained theoretically in terms of electron coherence using dynamic mean-field the- 
ory [37]. Using DFT, a model where spin-orbit interactions and orbital polarization 
are strong, but spin polarization is zero has been employed to yield a nonmagnetic 
configuration of 5-Pu [63]. This has received ancillary support by polarized neu- 
tron scattering measurements on single crystal PuCoGas, which shows the orbital 
moment dominates the magnetization [26]. Recent magnetic susceptibility measure- 
ments have shown that local magnetic moments in the order of 0.05 jllB/atom form 
in Pu as damage accumulates due to self-irradiation [38]. 



13 Conclusions 

The rare-earth metals have localized and atomic-like Af states across the series, 
leading to strong magnetic moments. The 5d transition metals are itinerant and 
band-like, behaving as typical metals with wide bands that strongly participate in 
bonding. The 5 / actinide metals exhibit a behavior similar as the 5d transition met- 
als in the early actinides Th, Pa, U, and Np, but then as the 4 / rare earths for the 
middle actinides Am, Cm, Bk, and Cf. This is directly attributed to a transition from 
itinerant to localized 5 / states that occurs at Pu. This transition can be examined 
through the N4.5 EELS spectra and the spin-orbit sum-rule analysis. The relative N4 
peak intensity reduces up to Am, then increases for Cm in EELS spectra, where the 
N4 and N5 peaks are only marginally broadened by multiplet structure. Across the 
actinide series, we see the light metals exhibit LS coupling while the middle metals 
Pu, Am, and Cm exhibit intermediate coupling. It is the transition from LS to inter- 
mediate coupling which reveals the transition from itinerant to localized 5/ states. 
This means that the EELS ^4,5 spectra and spin-orbit sum-rule analysis can be used 
as a measure of the degree of itinerancy in the actinide 5/ states. 

The EELS spectra and spin-orbit analysis clearly support that «/ ~ 5 and not 6 in 
Pu and that it falls near intermediate coupling curve. The 5 f 5 configuration with a 
spin-orbit interaction that adheres to intermediate coupling near the jj limit begs the 
question of why Pu is not magnetic. Recent advances in band theory have begun to 
address this question, such as DFT with zero spin polarization but strong spin-orbit 
and orbital polarization [63], DMFT calculations by Shim et al. [61] and Marianetti 
et al. [37], and LDA+U calculations by Shick et al. [60] 
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Americium falls precisely on intermediate coupling curve and has the largest 
spin-orbit interaction of all actinide metals. The results anchor the /-count for the 
adjacent actinides and clearly show the metal is non-magnetic because the 5 f 6 has a 
ground state 7 = 0. Finally, the results for Cm prove the spin-orbit sum rule works, 
showing the intermediate coupling curve bends back to LS curve. The experimental 
results [45] show that the angular-momentum coupling mechanism dictates large 
spin polarization and explain the magnetic stabilization of Cm observed by Heath- 
man et al. [23] 
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